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Abstract. The productions of real photons from quark gluon plasma and hot hadronic matter 
formed after the nucleus - nucleus collisions at ultra-relativistic energies are discussed. The effects 
of the spectral shift of the hadrons at finite temperature on the production of photons are investigated. 
On the basis of the present analysis it is shown that the photon spectra measured by WA98 collabora- 
tion in Pb + Pb collisions at CERN SPS energies can be explained by both QGP as well as hadronic 
initial states if the spectral shift of hadrons at finite temperature is taken into account. Several other 
works on the analysis of WA98 photon data have also been briefly discussed. 
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1. Introduction 

Calculations based on the QCD renormalization group predicts that strongly interacting 
systems at very high density and/or temperature are composed of weakly interacting quarks 
and gluons [1] due to asymptotic freedom and the Debye screening of colour charge. On 
the other hand at low temperature and density the quarks and gluons are confined within 
the hadrons. Therefore, a phase transition is expected to take place at an intermediate 
value of temperature and/or density. In fact such a transition, from hadronic matter to 
Quark Gluon Plasma (QGP) is actually observed in lattice QCD numerical simulations [2] 
at high temperature. One expects that ultra-relativistic heavy ion collisions (URHIC) at 
CERN/SPS, BNL/RHIC and CERN/LHC might create conditions conducive for the for- 
mation and study of QGP [3]. Among various signatures of QGP, photons and dileptons 
are known to be advantageous as these signals probe the entire volume of the plasma, 
with little interaction and thus, are better markers of the space-time history of the evolving 
matter [4]. 

The aim of the present work is to contrast the real photon emission rate from the follow- 
ing two nuclear collision scenarios: 



(1) A +A ^QGP^Mixed Phase^Hadronic Phase 
or 

(2) A + A ^Hadronic Phase, 
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by taking into account the finite temperature effects on the hadronic masses and decay 
widths. The main sources of photons from URHIC are: (a) the decay of mesons (mainly 
7T° and rf), (b) thermal source, either QGP and/or the hadronic reactions and hadronic de- 
cays in the thermal medium depending on which of the scenarios (1) or (2) realized after 
the collisions and (c) hard scattering of the partons embedded in the nucleons of the col- 
liding nuclei in the very early stage of the collisions. However, the transverse momentum 
(pr) spectra of single photons presented by WA98 collaboration [5] does not contain the 
contributions from the decays of mesons. Therefore, in the present article we will consider 
photons from (b) and (c) only. In the scenario (1) the pre-requisite for the QGP diagnostics 
(by measuring photon spectra) is to estimate the photon yield from hadronic sources and 
hard collisions of partons in the initial stage. 

In the next section we discuss the hard photon production from nucleus-nucleus colli- 
sions at SPS energies. Thermal photon emission from the QGP and hadronic matter are 
discussed in section 3. Hadronic properties at non-zero temperature is presented in section 
4. Section 5 is devoted to discuss the space time evolution and in section 6 we present our 
results and discussions. 



2. Hard QCD Photons 

The hard QCD photons are estimated using perturbative QCD as 

where Taa^) is the nuclear thickness at impact parameter b. Its value at b ~ 3.2 fm, 
corresponding to the most central event of WA98 experiment is <~ 220/fm 2 . a pp includes 
the pp cross-section for Compton and annihilation processes among the partons. At SPS 
energies one should include the effects of intrinsic kx distribution of partons [6] (due to 
finite size of the nucleons). This leads to substantial enhancement in the photon spectra [7]. 
In practice such an effect is implemented by multiplying each of the parton distribution 
functions appearing in the right hand side of the above equation by a Gaussian function of 
the type f(k T ) = exp[—k T /(h T )]/n(k T ) and integrating over d 2 kr- We use CTEQ5M 
partons [8] and (k T ) = 0.9 GeV 2 for evaluating the hard QCD photons. The energy at 
the centre of mass of Pb + Pb collisions at SPS is 17.3 GeV. Experimental data on hard 
photons does not exist at this energy. Therefore, the "data" at y/s = 17.3 GeV is obtained 
from the data at yfs = 19.4 GeV of the E704 collaboration [33] by using the scaling 
relation: Eda/d 3 p 7 \h a +h b -^x+i= f{%T = %Pt / \/s) / s 2 , for the hadronic process, 
h a + hf, — ► X + 7 [9]. However, such a scaling may be spoiled in perturbative QCD due 
to the momentum dependence of the strong coupling, a s and from the scaling violation of 
structure functions, resulting in faster decrease of the cross section than 1/s 2 . Therefore, 
the data at \fs = 17.3 GeV obtained by using the above scaling gives a conservative 
estimate of the prompt photon contributions. The photons from hard QCD processes have 
been used to reproduce the scaled p-p data of E704 collaboration. The higher order effects 
has been taken into account through a K-factor ~ 2. Now the question is, can we say 
that an enhanced production in A-A collisions compared to p-p will presumably mark the 
presence of a thermal source? Not necessarily, because in nucleus-nucleus collisions there 
may be enhancement in the high px part of the photon spectra due to various nuclear 
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effects, e.g., Cronin effects. However, the good news is that even if we take into account 
the transverse momentum broadening due to the finite size of the nucleon as well as due to 
the Cronin effects then we find that the theoretical yield is less than the experimental value 
(WA98 data) for 1.5 < p T (GeV)< 2.5 [10,11]. Clearly indicating the presence of thermal 
source. What is the nature of this thermal source which can reproduce the WA98 data? We 
will discuss this issue in the next section. 



3. Thermal photons from QGP and hadronic matter 

Since quarks are electrically charged their interactions at the thermal bath will produced 
photons we are looking for. The quark - anti-quark annihilation and QCD compton are 
the dominant processes for the production of photon from a thermalized system of quarks 
and gluons [12]. However, it has been shown [13] (see also [14,15]) that the two-loop 
contribution leading to bremsstrahlung and qq annihilation with scattering is of the same 
order as the lowest order processes. The total rate of emission (upto two loops) per unit 
four- volume at temperature T is given by 



dR 5 aa s 



In 



2.912 E 



; (Jr " j£) {ln2 + A] 



(2) 



where J T ~ 4.45 and J L ~ —4.26. The QCD coupling, 'g' is given by, 



= as = (33-2nf)H K T/T c ) 0) 

where rif is the number of quark flavours and k = 8 [16]. The photon emission rate in 
Eq. (2) is evaluated in the Hard Thermal Loop (HTL) approximation. However, the HTL 
approximation is not valid at SPS energies where an initial temperature about few hundred 
MeV may be realized. Because the HTL resummation is valid for g « 1 whereas the 
value of g obtained from Eq. (3) is <~ 2 at T <~ 200 MeV. At present it is unclear whether the 
rate in Eq. (2) is valid for such a large value of g or not. It would be rather difficult to make 
any firm conclusion from the results where Eq. (2) is used at SPS energies. Keeping this 
reservation in mind we will use Eq. (2) to evaluate the photon yield from QGP. However, 
we will see below that within the present framework the space time integrated photon yield 
from quark matter is less than that from hadronic matter due to the smaller life time of 
the QGP phase as a result of a moderate value of the initial temperature considered here. 
Therefore, the total thermal photon yield remains largely unaffected even in a scenario 
where QGP is formed in the initial state. 

The photon yield from the hadronic matter (HM) (see first of [12]), is evaluated from the 
reactions, n p — ► irj, nn — ► p 7, it it — > 77 7, 7r 77 — ► 7r 7 and the decays p — > 7r 7r 7 and 
uo — > ttj. We refer to Refs. [17] for the invariant amplitudes of these processes. Photon 
production due to the process tt p — > a\ — > 7r 7 is also considered here. 
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4. Hadronic properties at T ^ 

It has been emphasized that the hadronic properties will be modified due to its interactions 
with the particles in the thermal bath. As a consequence of the change in the properties 
of hadrons the static photon emission rates as well as the Equation of State (EOS) of the 
evolving matter will also change in a non-trivial way [18]. Broadly two kinds of medium 
modifications of hadrons are expected: (i) shift in the pole of the spectral function without 
any broadening and (ii) broadening of the spectral function but pole does not shift. In 
Ref. [18] the effects of spectral changes of hadrons on the electromagnetic probes is studied 
in detail. It was observed that the gauged linear and non-linear sigma models and the model 
with hidden local symmetry do not show any appreciable effect on photon emissions. In the 
Walecka model, the universal scaling hypothesis for the hadronic masses (except pseudo- 
scalar) has been seen to enhance photon emission. 

Both (i) and (ii) can reproduce the enhancement of the dilepton yield in the low invari- 
ant mass region, however, the photon yield is largely unaffected by (ii) since the spectral 
function is smeared out. Therefore, a simultaneous measurements of single photon and 
lepton pairs is important to shed light on the in-medium effects of hadrons. According to 
the universal scaling scenario [19] the in-medium quantities (denoted by *) at finite T is 
parametrized as 



mv 



fv 
fv 



U>0 



T 2 



(4) 



where V stands for vector mesons, fv is the coupling between the vector meson field and 
the electromagnetic current and ojq is the continuum threshold. Mass of the nucleon varies 
with temperature as in Eq. (4). It is to be noted that there is no definite reason to believe 
that all the in-medium dynamical quantities are dictated by a single exponent A. This is 
the simplest possible ansatz. The effective mass of a\ is estimated by the Weinberg's sum 
rules [20]. 

In the QHD model [21] the effective masses of nucleon, p and u mesons can be 
parametrized in the following forms: 



M* N = M N 



1 - 0.0264 



r(GeV) 
0.16 



■s.oi 



(5) 



771 „ 



m,„ = m u 



1-0.127 



1 - 0.0438 



T(GcV) 
0.16 

T(GeV) 
0.16 



5.24 



In (ii) the change of the width of the p meson with temperature is parametrized as, 

iy = r p /(i - t 2 /t c 2 ) 



(6) 



(7) 



and the mass remains constant to its vacuum value. The results for scenarios (i) and (ii) 
will be compared with a scenario (iii) where both the masses and widths of the hadrons 
remain fixed to their vacuum values. 
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Figure 2. Photon spectra at T = 180 MeV. Solid (dashed) line indicates result when 
hadronic masses vary according to Eq. (4) (fixed at vacuum values). Dotted line shows 
photon spectra from the decays of baryonic resonances. 
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5. Space Time Evolution 

It is assumed here that the produced matter reaches a state of thermodynamic equilibrium 
after a proper time ~ 1 fm/c [22]. In case of a deconfined matter is produced, it evolves 
in space and time till freeze-out undergoing a phase transition to hadronic matter in the 
process. The (3+1) dimensional hydrodynamic equations have been solved numerically 
by the relativistic version of the flux corrected transport algorithm [23], assuming boost 
invariance in the longitudinal direction and cylindrical symmetry in the transverse plane. 
The initial temperature Tj can be related to the multiplicity of the event, dN/ dy by virtue 
of the isentropic expansion as [24], 

dN 45C(3) 2 3 

-d-y- = ^2^* R ^ akT * n (8) 

where Ra is the initial radius of the system, Tj is the initial thermalization time and 
a-k = (7r 2 /90) gk', gk being the effective degeneracy for the phase k (QGP or hadronic 
matter). The bag model EOS is used for the QGP phase. gn(T), the statistical degeneracy 
of the hadronic phase, composed of ir, p, u>, r], a\ and nucleons is a temperature depen- 
dent quantity in this case and plays a crucial role in the EOS [18]. As a consequence the 
square of sound velocity, c~ 2 = \(T / gn){dgH / dT) + i] < 1/3, for the hadronic phase, 
indicating non-vanishing interactions among the constituents (see also [25]). The hydrody- 
namic equations have been solved with initial energy density, e(rj, r) [23], obtained from 
Ti through the EOS. We use the following relation for the initial velocity profile which has 
been successfully used to study transverse momentum spectra of hadrons [26,27], 

v ^^{-£) (9) 

Here we took 5 = 1 and the sensitivity of the results on u will be shown. It is observed 
that the results do not change substantially with reasonable variation of the parameter 5 for 
a given value of vq. 

The space-time integration can also be performed by taking the temperature profile from 
the transport model as [28,29], 

T(t) = {Ti - T co )e- T / T ° + (10) 

where t = 8 fm/c and T m = 120 MeV. 



6. Results and Discussions 



In Fig. (1) the effective masses of vector mesons is plotted as a function of temperature for 
universal scaling and QHD model calculations. In QHD the p and ui masses show different 
behaviour due to their different coupling strengths with the nucleons in the thermal bath. In 
Fig. (2), we show the photon emission rate from hadronic matter at a temperature, T = 180 
MeV. The solid and dashed lines correspond to the case (i) and (iii) respectively. The case 
with large collisional broadening shows no deviation from (iii). The increased photon 
yield at large energy (E) is caused by the enhancement in the Boltzmann factor due to the 
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Figure 3. The energy density e in the unit of T 4 for the equation of state used in 
the present work is plotted as function of temperature (T) in the unit of the critical 
temperature, T c . The filled circle denotes the lattice results [2]. 

reduction in meson (particularly p) masses. The thermal photon yield with hadronic mass 
variation due to Walecka model or Brown-Rho scaling [19] (A = 1/6 in Eq. (4)) will lie 
between the solid and dashed curves in Fig. (2). It is clearly observed that the contributions 
from the decays of baryonic resonances (7V(1520), iV(1535), iV(1440), A(1232), and 
A(1620)) are small (dotted line). The values of the decay widths, R — > N 7, where R 
and N denote the baryonic resonances and the nucleon respectively, are taken from the 
particle data book. 

We need the equation of state and the initial condition to solve the hydrodynamic equa- 
tions. The effects of the temperature dependent hadronic masses have been taken into 
account in the EOS through the effective statistical degeneracy [18]. In fig. 3 the temper- 
ature dependence obtained for different hadronic interactions is compared with the lattice 
QCD calculations [2]. The universal scaling scenario seems to reproduce the lattice data 
reasonably well. 

For central collisions of Pb nuclei at 158 AGeV at the CERN-SPS, we assume that 
QGP is produced at n=l fm/c which expands and undergoes a first order phase transi- 
tion to hadronic matter at T c =160 MeV. Taking dN/dy=700 and = <?qgp=37 for a 
two-flavour QGP, the initial temperature Tj comes out as 196 MeV. In a first order phase 
transition one has a mixed phase of coexisting QGP and hadronic matter which persists 
till the phase transition is over. Thereafter the hadronic matter expands, cools and freezes 
out at a temperature, Tf and radial velocity, v£ . The sum total of the photon yields from 
the QGP phase, the mixed phase and the hadronic phase constitutes the thermal yield [30]. 
The values of (Tf, v() should in principle be obtained from the analysis of hadronic spec- 
tra. In the present work the value of Tf is taken as 120 MeV which reproduces the hadron 
spectra [31] of NA49 collaboration [32]. The sensitivity of the results on the value of vq 
will be demonstrated below. 

In case of QGP formation the thermal photons contain contributions from quark matter 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 

P T (GeV) 



Figure 4. Total photon yield in Pb + Pb collisions at 158 A GeV at CERN-SPS. The 
theoretical calculations contain hard QCD and thermal photons. The system is formed 
in the QGP phase with initial temperature T» = 196 MeV. 

(QM = QGP + QGP part of mixed phase) and hadronic matter (HM = hadronic part of 
mixed phase + hadronic phase). In Fig. (4), results for the total photon emission is shown 
for three different values of the initial transverse velocity with medium effects as in case (i). 
All the three curves represent the sum of the thermal and the prompt photon contribution 
which includes possible finite k T effects of the parton distributions. The later, shown sep- 
arately by the dot-dashed line also explains the scaled pp data from E704 experiment [33]. 
We observe that the photon spectra for the initial velocity profile given by Eq. (9) with 
v = 0.3 explains the WA98 data reasonably well. It is found that a substantial fraction of 
the photons come from mixed and hadronic phases (hadronic masses vary with tempera- 
ture according to Eq. (4) with A = 1/2) . The contributions from the QGP phase is small 
because of the small life time of the QGP (~ 1 fm/c). Therefore, the results shown in this 
figure is largely independent on the uncertainties (mentioned above) involved in the photon 
emission rates from QGP. 

The above statement together with the uncertainty of the critical temperature T c [2] 
poses the following question: Is the existence of the QGP phase essential to reproduce 
the WA98 data? To study this problem, we have considered two possibilities: (iii) pure 
hadronic model without medium-modifications, and (i) pure hadronic model with scaling 
hypothesis according to Eq.(4) (for A = 1/2). In the former case, Ti is found to be <~ 250 
MeV for Tj = 1 fm/c and dN/dy — 700, which appears to be too high for the hadrons 
to survive. Therefore we exclude this possibility here. On the other hand, the second 
case with an assumption of Ti = T c (which is just for simplicity) leads to Ti <~ 205 
MeV, at Ti = 1 fm/c, which is not unrealistic. In this case, the hadronic system expands 
and cools and ultimately freezes out at T/=120 MeV. The masses of the vector mesons 
increase with reduction in temperature (due to expansion) according to Eq.(4). The results 
of this scenario for three values of the initial radial velocity including the prompt photon 
contribution are shown in Fig. (5). The experimental data are well reproduced for vanishing 
initial transverse velocity also. Therefore, the results shown in Fig. (4) and (5) indicate that 
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a simple hadronic model is inadequate. Either substantial modifications of hadrons in the 
thermal bath or the formation of QGP in the initial stages is necessary to reproduce the 
data. It is rather difficult to distinguish between the two at present. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 

P T (GeV) 



Figure 5. Same as Fig. (4) with hadronic initial state at T; = 200 MeV. 



In Fig. (6) we compare the results of hydrodynamic and transport calculations. Within 
the framework of the transport model the data is well reproduced when the hadronic masses 
are allowed to vary (scenario i) according to the Eq. 4 (A = l/2,long dashed line). As 
mentioned before the photon spectra is insensitive to the change in width (scenario ii). 
In the scenario (ii) the experimentally observed "excess" photon in the region 1.5 < pr 
(GeV) < 2.5 (dotted line) is not reproduced. The change in the width of the vector mesons 
(p in particular) has very little effects on the pt spectra of photons due to the following 
reasons. The density of an unstable particle in a thermal bath can be written as [34], 

</V " e-^/ T P(s) (11) 



d 3 kd 3 xds (2tt) 

where g is the statistical degeneracy of the particle and P(s) is the spectral function, 

1 Imll 

7T (s - m 2 - Ren) 2 + (ImlT) 



ns) = - - ^ 2 p „ m2 , , T „ m2 d2) 



Imll (Re IT) is the imaginary (real) part of the (trace of) p self energy. Eqs. (11) and 
(12) indicate that the density of particles in a thermal bath is given by the Boltzmann 
distribution weighted by the Breit-Wigner function, which gets maximum weight from the 
value of s = m 2 + Re IT, the contribution from either side of the maximum being averaged 
out. Therefore, the results become sensitive to the effective mass, s = m 2 * = m 2 + RelT 
and not to the width of the spectral distribution. The dash-dotted line indicates results with 
vacuum masses and widths. In case of transport model calculations (dashed line) there 
is excess photons at the low px region compared to the hydrodynamic model (solid line) 
due to the following reason. We find that the variation of temperature with time (cooling 
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law) in Eq. 10 is slower than the one obtained by solving hydrodynamic equations. As a 
consequence the thermal system has a longer life time than the former case, allowing the 
system to emit photons for a longer time. In case of hydrodynamics this is compensated 
by the transverse kick experienced by the photon at large pt due to radial velocity of the 
expanding matter. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 

P T (GeV) 

Figure 6. Same as Fig. (5). Results for hydrodynamic and transport (indicated by 
'trans') model calculations are compared. 

It is widely believed that the photon spectra can be very useful for the estimation of 
the initial temperature. In table I We show the initial temperature obtained by several au- 
thors [35,14,36-38] by analyzing the WA98 photon spectra. In [37] the thermal photon 
spectra was parametrized as EdN/d 3 p = V±F{E, T e ff,v ), similar to the parametrization 
one uses to study the thermal hadronic spectra. A value of T e ff = 170 MeV and v = 0.3 
can reproduce the WA98 data if the prompt photon contribution is normalized to reproduce 
the data forp^ > 2 GeV. As T e f f is the average value of the temperature (Ti > T e f / > Tf) 
in this analysis, one may expect an initial temperature ~ 200 MeV. In [38] the hard pho- 
ton contribution has been used to normalized the scaled p-p data at ^fs = 19 GeV and the 
thermal photon spectra was evaluated for both QGP and hadronic initial state assuming a 
non-zero radial velocity of the QGP/hadronic fluid. Steffen and Thoma [14] have demon- 
strated that the value of the statistical degeneracy (gh) in the hadronic phase (and hence 
the life time of the mixed phase) is crucial for the description of WA98 data. An initial 
temperature ~ 220 MeV with T c = 170 MeV and gh = 8 can reproduce the data in this 
case. In [35] and [14] the static photon emission rate from QGP is similar, then why the 
value of Ti is so different? Apparently the main reason is the difference in the life time of 
the mixed phase between the two cases. In [14] the life time of the mixed phase is large 
(small hadronic degeneracy) compared to that of ref. [35], hence allowing the mixed phase 
to emit photon for a longer time interval in the former case. In [14] the contributions from 
the QGP phase due to smaller initial temperature (compared to [35]) is compensated by 
a larger contributions from the mixed phase. Huovinen et. al. [36] has studied in detail 
the effects of various EOS, boost invariant and non-invariant hydrodynamic flow on the 
photon spectra. Thermal photon with values of 7* <~ 213 — 255 MeV (depending on the 
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EOS and evolution scenario) + hard QCD photon can reproduce the data well in their case. 
We would like to mention here that the (thermal) photon emission rate used in [36] is full 
order of a s results obtained in [39]. 



[35] 


[37] 


[14] 


[38] 


[36] 


Present 


Ti 


Ti 


Ti 




T 


T 


335 


210 


220 


200-230 


213-255 


200 



Table I: The initial temperature obtained in the present analysis is compared with the values 
obtained by other authors [35,37,14,38,36]. 

In order to reproduce the WA98 photon data either a substantial reduction in vector 
meson masses or the formation of QGP in the initial stage with Tj <~ 200 MeV is necessary. 
A simple hadronic model is appears to be ruled out by the experimental data. 

In spite of the encouraging situation mentioned above, a firm conclusion about the for- 
mation of the QGP at SPS necessitates a closer look at some pertinent but unsettled issues. 
A prerequisite for the detection of the QGP by studying the photon spectra is to subtract 
contribution from the initial hard processes. Therefore, it is extremely important to know 
quantitatively the contribution from the hard processes. Again, the assumption of com- 
plete thermodynamic equilibrium for quarks and gluons may not be entirely realistic for 
SPS energies; lack of chemical equilibrium of quarks will further reduce the thermal yield 
from QGP. We have assumed Tj = 1 fm/c at SPS energies, which may be considered as the 
lower limit of this quantity, because the transit time (the time taken by the nuclei to pass 
through each other in the CM system) is <~ 1 fm/c at SPS energies and the thermal system 
is assumed to be formed after this time has elapsed. In the present work, when QGP initial 
state is considered, we have assumed a first order phase transition with bag model EOS 
for the QGP for its simplicity, although it is not in complete agreement with the lattice 
QCD simulations [2]. However, it is difficult to distinguish among different EOS with the 
current resolution of the photon data [36]. As mentioned before, there are uncertainties 
in the value of T c [2], a value of T c <~ 200 MeV may be considered as an upper limit. 
Moreover, the photon emission rate from QGP given by Eq. (2), evaluated in Refs. [12,13] 
by resumming the hard thermal loops is strictly valid for g < < 1 whereas the value of g 
obtained from Eq. (3) is <~ 2 at T <~ 200 MeV. At present it is not clear whether the rate in 
Eq. (2) is valid for such a large value of g or not. New method is required to evaluate the 
photon emission rate from QCD Plasma, HTL approximation is not valid at SPS energies 
and it may not be valid even at LHC energies. Evaluation of the spectral function from 
lattice QCD [40] will be very useful in this context. 
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